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ABSTRACT 

The isophotal wavelengths, flux densities, and AB magnitudes for Vega 
(a Lyr) are presented for the Mauna Kea Observatories near-infrared filter set. 
We show that the near-infrared absolute calibration for Vega determined by Co- 
hen et al. and Megessier are consistent within the uncertainties, so that either 
absolute calibration may be used. 

Subject headings: instrumentation: photometers — techniques: photometric 



1. INTRODUCTION 

Simons & Tokunaga (2002) and Tokunaga, Simons, & Vacca (2002) defined a 1-5 jum 
filter set that is designed to maximize sensitivity while minimizing the effects of atmospheric 
absorption on reducing the signal-to-noise ratio. This filter set is intended to provide good 
transformation between observatories located at altitudes of 2-4 km. Filter production runs 
have been organized to produce these filters, which are in use at more than 30 institutions. 
This filter set provides greater transmission than those advocated by Young, Milone, & Stagg 
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(1994), and was designed to provide nearly ideal photometric accuracy. To distinguish this 
filter set from others, we refer to it as the Mauna Kea Observatories near-infrared (MKO- 
NIR) filter set. 

In this paper, we present the isophotal wavelengths, flux densities, and AB magnitudes 
for Vega for the MKO-NIR filter set. We compare the absolute calibration advocated by 
Cohen et al. (1992) to that of Megessier (1995), and we show that there is no significant 
difference in the flux densities for Vega derived by these authors. 



2. ISOPHOTAL WAVELENGTHS AND ZERO MAGNITUDE FLUX 

DENSITIES 



The number of photo-electrons detected per second from a source with an intrinsic 
spectral energy distribution F\{X) is given by 



Here T(A) is the atmospheric transmission, Q(X) is the product of the throughput of the 
telescope, instrument, and quantum efficiency of the detector, R(X) is the filter response 
function, and A te \ is the telescope collecting area. The system response S(X) is equal to the 
relative spectral response (RSR) defined by Cohen et al. (2003). 

If F\(X) and S(X) are both continuous and S(X) is nonnegative over the wavelength 
interval, then from equation (2) and the mean value theorem for integration there exists a 
A iso such that 



2.1. The Definition of Isophotal Wavelength 




(1) 



(2) 



where S(X) is the total system response given by 



S(X) = T(X)Q(X)R(X)A tel 



(3) 





Rearranging this, we obtain 



(5) 
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where A iso denotes the "isophotal wavelength" and (Fx) denotes the mean value of the in- 
trinsic flux above the atmosphere (in units of W m -2 /im" 1 ) over the wavelength interval of 
the filter. Thus Ai SO is the wavelength at which the monochromatic flux F\(\± so ) equals the 
mean flux in the passband. Hence Ai SO and -Fa^so) are the wavelength and monochromatic 
flux density, respectively, that best represent a broadband (heterochromatic) measurement. 
In addition we choose to use isophotal wavelengths for consistency with the extensive series 
of papers on infrared calibration by Cohen and collaborators. 



In a similar fashion, 



F v {y*>) - (F v ) - j s{v)/vdv , (6) 



where u- mo denotes the "isophotal frequency" and (F u ) denotes the mean value of the intrinsic 
flux above the atmosphere (in units W m~ 2 Hz -1 ) over the frequency interval of the filter. 

We show in Table 1 the MKO-NIR filter isophotal wavelengths for Vega. These were 
calculated from equation (5) using an atmospheric model of Vega computed by R. Kurucz 1 
with the parameters T cS = 9550 K, log(g) = 3.95, v tuib = 2 km s" 1 , v rot = 25 km s _1 , and 
[Fe/H] = —0.5. These parameters are the same as those adopted by Bohlin & Gilliland 
(2004). The model has a resolving power of 10 5 and has been scaled to the absolute flux 
level of 3.46 x 10~ 9 erg cm -2 s" 1 A" 1 at 5556 A as determined by Megessier (1995). We 
used ATRAN (Lord 1992) to calculate the atmospheric transmission T(A) for a range of 
precipitable water values between and 4 mm, an airmass of 1.0, and the altitude of Mauna 
Kea. The measured filter response curves 2 given by Tokunaga et al. (2002) were used for 
-R(A) (see their Fig. 1). We assumed that the throughput term Q(X) was a constant over the 
wavelength integrals. Thus our calculations are precisely correct only for constant detector 
responsivity and instrumental throughput. However Stephens & Leggett (2004) find that 
the variations in the detector responsivity and instrumental throughput leads to photometric 
variations of <0.01 mag. This corresponds to variations of A iso of <0.3%. 

In accordance with the design goals of the MKO-NIR filters, variations in A iso as a 
function of precipitable water vapor were found to be very small (<1%) over the range of 
water vapor values typically encountered on Mauna Kea. For the range of to 4 mm of 
precipitable water vapor, we found variations AAj SO of 0.000 /im (J), 0.001 /im (H), 0.012 
/xm (K'), 0.013 /im (K s ), 0.001 /im (K), -0.014 /mi (L'), and 0.001 /im (M'). 

The methods employed here, in particular using the number of photons detected and 



1 http://kurucz.harvard.edu/stars/VEGA 

2 http://irtfwcb. ifa.hawaii.edu/Facility/nsfcam/hist/newfilters. html 
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the formulation in equation (5), which includes the wavelength terms in the integrals, are 
the same as those used by Cohen et al. (1992) and subsequent papers by Cohen and his 
collaborators. For clarity we have explicitly presented the equations for calculating the 
isophotal wavelengths for photon counting detectors. See Bessell et al. (1998) for discussion 
of how photometric results differ between energy measuring detectors and photon-counting 
detectors. 

Real spectra do not necessarily satisfy the requirements of the mean value theorem 
for integration, as they exhibit discontinuities. Although the mean value of the intrinsic 
flux is well-defined, the determination of the isophotal wavelength becomes problematic 
because real spectra contain absorption lines and hence the definition can yield multiple 
solutions. In addition, A iso is not an easily measured observational quantity because it 
depends on knowledge of the intrinsic source flux distribution F\(\), which is exactly what 
one is attempting to determine with broadband photometry. Nevertheless for the reasons 
stated above we use isophotal wavelengths in this paper. Other definitions of the filter 
wavelength are briefly discussed in the Appendix for completeness. 

Any definition of the filter wavelength suffers from the limitation that the spectral energy 
distribution of the object being observed is likely to be different from that of Vega. Since the 
isophotal wavelength is different for objects that have different spectral energy distributions, 
a correction factor is required to obtain the monochromatic magnitude at the same isophotal 
wavelength as Vega. Hanner et al. (1984) discuss this problem in detail for observations of 
comets. 



2.2. Absolute Flux Densities for Vega 

The question of the near-infrared absolute flux densities for Vega above the atmosphere 
has been discussed by Cohen et al. (1992) and Megessier (1995). Cohen et al. determined 
their absolute calibration from a model atmosphere for Vega multiplied by the atmospheric 
transmission and instrument response (filters, throughput, and filter response). They did not 
use absolute calibration measurements by Blackwell et al. (1983) and Selby et al. (1983) be- 
cause Blackwell et al. (1990) concluded that atmospheric models of Vega offered higher pre- 
cision than the observationally determined absolute calibration in the near-infrared. Bessell 
et al. (1998) also concluded that model atmospheres are more reliable than the near-infrared 
absolute calibration measurements. However Megessier (1995) argued that the models were 
not reliable and that the near-infrared absolute calibration of Vega should be based on mea- 
surements that are independent of atmospheric models. Based on four model-independent 
measurements, Megessier (1995) determined an averaged absolute flux density for Vega. 



- 5 - 



We directly compared the absolute calibration of Vega determined by Cohen et al. 
(1992) and Megessier (1995). We first used a third-order polynomial to fit the logarithm of 
the flux density for Vega as a function of the logarithm of the wavelength from Megessier. 
The use of logarithms gives a nearly linear relationship. The wavelengths and flux densities 
for 0.0 mag were taken from Table 4 in Megessier. We assumed a near-infrared magnitude 
of 0.02 mag for Vega to be consistent with the visible magnitude assumed by Megessier. We 
then interpolated to the wavelengths cited by Cohen et al. in their Table 1 to make direct 
comparisons to the Megessier results. 

The difference in the logarithm between the Megessier and Cohen et al. absolute flux 
densities for Vega are shown in Figure 1. The mean of the difference is —0.0022 ± 0.0031. 
Thus the absolute flux densities of Vega as given by Megessier and Cohen et al. are indis- 
tinguishable. The uncertainty of the Megessier and Cohen et al. flux densities for Vega are 
about 2% and 1.45%, respectively. 

While the results discussed above show the consistency of the independent methods of 
Cohen et al. and Megessier, we note the following caveats: 

1. Both Cohen et al. (1992) and Megessier (1995) rely on the absolute calibration of Vega 
at 0.5556 fj,m. Cohen et al. adopt a flux density of 3.44 x 10" 8 W m~ 2 //m -1 , while 
Megessier adopts 3.46 x 10~ 8 W m~ 2 /im" 1 . Thus Megessier's flux density for Vega is 
0.6% higher at V. 

2. Recent work by Gulliver et al. (1994) and Peterson et al. (2004) indicate that Vega is 
pole-on and a fast rotator. Thus standard model atmospheres are not appropriate for 
Vega as discussed by Bohlin & Gilliland (2004). Nonetheless we have already shown 
the agreement between the Cohen et al. values and the model-independent results of 
Megessier. In addition, Price et al. (2004) show that the Midcourse Space Experiment 
(MSX) absolute calibration experiment is in agreement with the Cohen et al. (1992) 
values for Vega, to within the experimental errors of 1%. 

We show in Table 1 the flux densities for Vega for the isophotal wavelengths of the 
MKO-NIR filters. We have adopted the 1-5 /im flux densities for Vega as presented by 
Cohen et al. (1992) in their Table 1. We first computed the flux densities for Vega assuming 
a precipitable water vapor value of 2 mm at an airmass of 1.0, in addition to the parameters 
for Vega discussed in Section 2.1. Since the model flux density for Vega is 1.9% lower than 
that used by Cohen et al. (1992), we increased our calculated values by 1.9%, and this is 
shown in Table 1. Bohlin & Gilliland (2004) also found that the flux density of Vega in the 
infrared was about 2% lower than that presented by Cohen et al. (1992). 
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For the V isophotal wavelength and flux density calculations, we used the absolute 
spectrophotometry for Vega given by Bohlin & Gilliland (2004). The V filter profile used 
was that of Landolt (1992), obtained from Cohen et al. (2003; see the electronic version of 
the paper). The values we obtained for this V filter are shown in Table 1. 

It is evident from the above discussion that there is no consistent published atmospheric 
model for Vega at both visible and infrared wavelengths. This is primarily because Vega is 
nearly pole-on and there is a range of temperature from the hotter pole regions to the cooler 
equatorial regions. A single temperature model for Vega is therefore not realistic. The 
values of isophotal wavelengths and flux densities shown in Table 1 represent a best estimate 
based on absolute calibrations using blackbody sources, observations of standard stars, and 
atmospheric models. 

We note that Megessier (1995) found that the observed fluxes of Vega were about 2% 
higher than the atmospheric models. This problem was attributed to a possible near-infrared 
excess of Vega. However Leggett et al. (1986) found no infrared excess from Vega compared 
to other AO stars. This suggests that atmospheric models for Vega at near-infrared wave- 
lengths are in error, possibly because Vega is observed pole-on. 

2.3. Comment on the Definition of Zero Magnitude 

Infrared photometric systems at 1-5 /im are usually defined as being based on the John- 
son system or in a system in which the magnitude of Vega is taken to be 0.0. Examples of 
the former include systems at the Univ. of Arizona (Campins et al. 1985), ESO (Wamsteker 
1981), SAAO (Carter 1990), and AAO (Allen & Cragg 1983). In these systems, the magni- 
tude of Vega is defined as 0.02 or 0.03 mag. Examples of the latter include the systems CIT 
(Elias et al. 1982) and the Las Campanas Observatory (Persson et al. 1998). The UKIRT 
photometric system (Hawarden et al. 2001; Leggett et al. 2003) is based on the Elias et al. 
(1982) standard stars, so it follows the convention that the magnitude of Vega is 0.0 mag. 
Cohen et al. (1992) adopted a magnitude of 0.0 mag at infrared wavelengths, so that the flux 
density of Vega defines the flux density for 0.0 mag. This is continued in subsequent papers, 
and in Cohen et al. (2003) the nonzero magnitude of Vega at optical wavelengths is taken 
into account. Thus when applying the results in Table 1, one must take into consideration 
which photometric system is used. 
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2.4. AB Magnitudes 

The monochromatic AB magnitudes were defined by Oke & Gunn (1983) as 

AB = -2.51og(f„) -48.60 (7) 

where f v is in units of ergs cm -2 s -1 Hz -1 (see also Fukugita et al. 1996). The constant 
is set so that AB is equal to the V magnitude for a source with a flat spectral energy 
distribution. We adopt the Vega flux densities recommended by Bohlin & Gilliland (2004) 
(their alpha_lyr_stis_002.fits file). The visible flux values are tied to a flux density of 3.46 x 
10 -8 W m -2 yum -1 at 0.5556 fim following Megessier (1995). The isophotal wavelength at 
V is 5546 A 3 , and the isophotal flux density is 3.63 x 10 -20 ergs s _1 cm" 2 Hz -1 (Table 1). 
A V magnitude of 0.026 from Bohlin & Gilliland (2004) is assumed. Then 

AB = -2.5 log(f^) - 48.574 (8) 

or for F„ expressed in units of Jy, we have 

AB = -2.5 log(F„) + 8.926 . (9) 

The AB magnitudes for Vega were calculated from equation (9) and are shown in the AB 
magnitude column of Table 1. A great advantage of AB magnitudes is that the conversion 
to physical units at all wavelengths can be obtained with a single equation: 

F v = 3720 10~ OAAB . (10) 

The constant in equation 8 differs from that of Fukugita et al. (1996) because for Vega 
we assumed a different flux density value at V and adopted a different visual magnitude. 
However, it is within the uncertainty of the absolute calibration of 2% for Vega stated by 
Oke & Gunn (1983). 



3. SUMMARY 

1. Isophotal wavelengths, flux densities, and AB magnitudes for Vega are derived for 
MKO-NIR filter set. 

2. The absolute calibrations by Cohen et al. (1992) and Megessier (1995) are shown to be 
identical within the uncertainties. We adopt the 1-5 /xm absolute calibration of Cohen 
et al. to be consistent with the subsequent papers by Cohen and his colleagues. 



3 Should be 5450 A, sec erratum. 



- 8- 



3. The \^-band isophotal wavelength and flux density is given for completeness using the 
renormalized Vega model and the absolute calibration adopted by Bohlin & Gilliland 
(2004). The constant in the AB magnitude definition was determined from the recent 
Hubble Space Telescope (HST) measurements of the V magnitude and flux density of 
Vega and differs from that defined by Oke & Gunn (1983). 

4. There is no self-consistent atmospheric model for Vega at visible and infrared wave- 
lengths. Improvements to Table 1 can be expected with models that take into account 
that Vega is observed pole-on, and also with observations of a grid of AO stars (includ- 
ing Sirius) to eliminate the dependence on currently unreliable atmospheric models for 
Vega in the infrared. 



We thank M. Cohen, D. Peterson, and T. Nagata for useful discussions and Steve 
Lord for making ATRAN available to us. A.T.T. was supported by NASA Cooperative 
Agreement number NCC 5-538. This research has made use of NASA's Astrophysics Data 
System Bibliographic Services. 



4. APPENDIX A 

Other definitions for effective wavelengths are briefly discussed here. A more detailed 
discussion may be found in Golay (1974). In a manner similar to that used in equation (4), 
we can define the "effective wavelength" as 

A eff J F(X)S(X)d\ = J \F\(\)S(\)d\ , (Al) 

and hence 

Ks ~ jF x (X)S(X)dX ■ (A2) 

This is the wavelength at which the flux distribution in energy units, integrated over 
the passband and then converted to photons, equals the flux distribution in photon units 
integrated over the passband. Alternatively, F\(\)S(\) can be thought of as a probability 
distribution for the detection of energy from the source, and therefore, A e fr is the mean 
wavelength of the passband as weighted by the energy distribution of the source over the 
band. In a similar manner, we can determine a mean wavelength for the passband as weighted 
by the photon distribution of the source over the band. In this case, 
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_ J X 2 Fx(X)S(X)d\ 
Ks ~ j XF x (X)S(X)dX • [A6) 

Here, XF X (X)S(X) can be thought of as the probability distribution for detecting a 
photon from the source. Note that both effective wavelength definitions depend on the 
spectral energy distribution of the source. 

We can define a source-independent wavelength as follows: 

J P x {X)S{X)dX = J F x {X)S{X)/hvdX , (A4) 
where P\{X) = F\(\)/hv is the photon flux from the object. Thus, 

(P A ) J S{X)dX = ^ J XS(X)dX . (A5) 
Setting (Pa) = X (F x )/hc yields 

_ JXS(X)dX 

Ao " JS(X)dX ' (A6) 
which is the "mean wavelength" of the system. 

We can express equation (1) in terms of frequency to derive another filter wavelength. 
From equation (1) we have 

N p =y c J ^Fx(X)S(X)dX = i J F v {y)S{y)dvjv . (A7) 

In a manner similar to the derivation given above for equation (14), 

Va) J XS(X)dX = (F„) J ^-du (A8) 

= (F v ) j ^dX , (A9) 
where the last equation results from the fact that dvjv = dX/X. If we set 

(F v ) = (F x }Xl ivot /c , (A10) 

we obtain 

/ / XS(X)dX 

A pivot - y j s(x)dx/x ' (A11) 

which is known as the "pivot wavelength" of the system. The pivot wavelength provides an 
exact relation between F v and Pa given by equation (A10). This definition is used in the 
HST Synphot Users Guide (Bushouse & Simon 1998; see also Koorneef et al. 1986). 

The wavelengths based on the different definitions are shown in Table 2. 
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Fig. 1. — Difference of the logarithm of the Cohen et al. (1992) and Megessier (1995) flux 
densities for Vega. The Megessier flux density was fitted with a third order polynomial and 
subtracted from the Cohen et al. values (see text). 
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5. ERRATUM - to be submitted 

The isophotal wavelength for the V filter given in Section 2.4 and Table 1 should be 
5450 A. 

We did not show the isophotal frequency in Table 1 that should be used with the 
flux density in frequency units (Jy) and the AB magnitudes derived from it. Therefore 
we present in Table 3 the isophotal frequencies. Note that the isophotal wavelength and 
isophotal frequency must be derived using equations (5) and (6) since \ iso ^ c/u iso . In 
addition, F„ and AB magnitudes should always be plotted with v iso . 
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Table 1. Isophotal wavelength, flux densities, and AB magnitudes for Vega 



Filter 


Aiso 


F A 


F„ 


AB magnitude 




(H 


(W m~ 2 /iin -1 ) 


(Jy) 




V 


0.5446 


3.68E-08 


3630 


0.026 


J 


1.250 


3.01E-09 


1560 


0.943 


H 


1.644 


1.18E-09 


1040 


1.38 


K' 


2.121 


4.57E-10 


686 


1.84 


Ks 


2.149 


4.35E-10 


670 


1.86 


K 


2.198 


4.00E-10 


645 


1.90 


V 


3.754 


5.31E-11 


249 


2.94 


W 


4.702 


2.22E-11 


163 


3.40 



Note. — We assume the Landolt V filter profile (see text) 
and 2mm of precipitable water. 



Table 2. Isophotal, Effective, Mean, and Pivot Wavelengths for the MKO-NIR filters 



Filter 


Aiso 


Aeflf 


Aeflf 


Ao 


Apivot 




(H 


(/iin) 




(urn) 


(/iin) 


J 


1.250 


1.241 


1.243 


1.248 


1.247 


H 


1.644 


1.615 


1.619 


1.630 


1.628 


K' 


2.121 


2.106 


2.111 


2.123 


2.121 


Ks 


2.149 


2.138 


2.141 


2.151 


2.150 


K 


2.198 


2.186 


2.190 


2.202 


2.200 


L' 


3.754 


3.717 


3.727 


3.757 


3.752 


M' 


4.702 


4.680 


4.681 


4.684 


4.684 
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Table 3. Isophotal frequency, flux densities, and AB magnitudes for Vega 



filter 






AB mag 


name 


(xlO 14 Hz) 


(Jy) 


(mag) 


V 


5.490 


3630 


0.026 


J 


2.394 


1560 


0.941 


H 


1.802 


1040 


1.38 


K' 


1.413 


686 


1.84 


Ks 


1.395 


670 


1.86 


K 


1.364 


645 


1.90 


V 


0.7982 


249 


2.93 


W 


0.6350 


163 


3.40 



0.02 - 



0.00 - 



-0.02 - 



-0.04 - 



-0.06 



wavelength (/xm) 
1.5 2.0 3.0 




0.0 0.1 



0.2 0.3 0.4 0.5 
log wovelength (//m) 



0.6 0. 



